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The  work  described  in  this  publication  was  performed  by  the 
Mathematical  Analysis  Research  Corporation  (MARC)  under  contract  to 
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Institute  of  Technology.  This  activity  is  sponsored  by  the  Jet 
Propulsion  Laboratory  under  contract  NAS7-918,  RE182,  A187  with  the 
National  Aeronautics  and  Space  Administration,  for  the  United  States 
Army  Intelligence  Center  and  School. 

This  specific  work  was  performed  in  accordance  with  the  FY-87 
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I. 


INTRODUCTION 


The  purpose  of  this  report  is  to  analyze  the  various  fix  algorithms  used 
in  determining  estimates  of  emitter  location  (fixes)  based  on  varying  numbers 
of  Lines  of  Bearing  (LOBs).  Since  sensors  do  not  produce  infallible  readings, 
observations  (LOBs  on  an  emitter)  are  subject  to  measurement  error  and  the 
actual  estimated  location  will  differ  from  the  true  location.  The  difference 
between  the  expected  value  of  the  estimates  and  the  true  emitter  location  is 
called  bias1 .  This  analysis  focuses  on  how  the  bias  of  the  different 
algorithms  is  affected  by  increasing  or  decreasing  sample  size  (number  of 
LOBs),  varying  sensor  error,  and  changing  the  relative  locations  of  emitters 
and  sensors. 

From  observations,  it  is  possible  to  calculate  an  estimate  of  the  true 
location  of  the  emitter  by  utilizing  different  fix  algorithms.  Four 
2-dimensional  algorithms  were  chosen  for  study  in  this  report:  the 
Perpendicular,  the  Weighted  Perpendicular ,  the  Sine  of  Error  Minimization,  and 
the  Minimization  of  Angular  Error  methods.  These  algorithms  were  chosen 
because  they  represent  variations  on  classical  fix  theory  and  because  they 
relate  to  fix  methods  of  actual  systems.  For  instance,  the  Perpendicular 
method  is  the  'classical'  method  introduced  by  Stansfield;  the  Weighted 
Perpendicular  method  is  similar  to  the  fix  algorithm  in  Guardrail  V;  the  Sine 
of  Error  Minimization  method  is  comparable  to  that  used  in  FFIX;  and  the 
Minimization  of  Angular  Error  method  is  the  theoretically  'ideal'  approach 
according  to  some  sources.  Although  they  are  indeed  variations,  the  four 
methods  presented  here  are  not  altogether  different  from  each  other. 

In  general,  each  of  the  four  fix  methods  mentioned  above  is  biased. 
General  expressions  (located  in  the  appendix)  for  first-order  bias*  have  been 
derived  for  each  of  the  methods.  In  their  general  form,  these  expressions  can 
be  used  for  calculations  but  are  too  complicated  to  interpret  qualitatively. 
So,  the  qualitative  analysis  of  bias  has  been  divided  into  Intermediate  steps. 
Since  the  bias  of  the  methods  changes  as  a  function  of  the  number  of  LOBs, 
bias  has  first  been  analyzed  at  the  extreme  bounds  of  sample  size:  2  LOBs  and 
infinite  LOBs.  Then,  special  cases  have  been  examined  between  these  bounds  to 
obtain  a  better  Intuitive  grasp  for  the  behavior  of  bias. 

Conclusions  drawn  from  these  expressions  are  supported  by  simulation 
results  and  will  be  illustrated  by  diagrams.  It  is  the  goal  of  this  report  to 
examine  certain  pertinent  questions  regarding  bias.  Which  situations  or  which 
fix  methods  produce  the  most  bias?  How  sensitive  is  bias  to  sample  size  and 
other  factors?  Answers  to  these  and  similar  questions  are  discussed  in  the 
following  text. 


1  Formally,  let  (x,y)  be  the  estimated  location,  based  on  a  particular 
algorithm,  and  (x0,y0)  be  the  true  location.  Then,  bias  is  defined  to  be 

A  A 

the  vector  from  (x0,y0)  to  the  expected  value  of  (x,y). 

*  Bias  has  been  expanded  using  a  multivariate  Taylor  Series  expansion.  The 
higher  order  terms,  which  are  generally  insignificant,  have  been  dropped.  See 
the  appendix  for  more  detail. 


II.  CONCLUSION 


& 


The  particular  properties  for  the  various  methods  in  special  cases  are 
discussed  in  depth  in  section  III.  Here,  however,  broad  results  concerning 
bias  are  summarized.  First,  results  pertaining  to  the  behavior  of  bias  in 
general  are  outlined.  Then,  the  way  in  which  bias  behaves  as  a  function  of 
sample  size  has  been  carefully  described.  For  detailed  descriptions  of  the 
various  methods  and  how  they  operate,  see  the  appendix. 

For  the  purposes  of  this  report,  bias  has  been  defined  to  be  a  vector. 

As  such,  it  has  a  magnitude  and  a  direction.  Both  the  magnitude  and  the 
direction  are  affected  by  sample  size  and  the  configuration  of  the  sensors. 

In  addition,  however,  the  magnitude  is  also  affected  by  the  standard  deviation 
of  the  angular  errors. 

The  standard  deviation  (a)  associated  with  the  sensor  measurement  error 
will  affect  the  size  of  bias,  regardless  of  the  method,  the  configuration,  or 
the  sample  size.  The  magnitude  of  bias  is  proportional  to  o*.  In  Figures 
1,2, and  3.  the  large  region  surrounding  the  emitter  is  an  Indication  of  the 
spread  of  the  estimates  while  the  smaller  region  in  contact  with  the  emitter 
represents  bias.  Notice  that  when  a  has  been  doubled,  the  region  of  bias  has 
become  four  times  larger.  In  Figure  3.  a  has  been  tripled  and,  subsequently, 
the  region  of  bias  has  increased  by  a  factor  of  nine.  Thus,  the  extent  to 
which  bias  can  be  a  problem  clearly  depends  heavily  on  the  size  of  o. 

Another  major  factor  in  determining  the  magnitude  of  bias  is  the 
configuration  of  the  sensors.  The  simplest  situation  which  will  produce  large 
bias  is  any  configuration  In  which  the  angles  between  the  true  lines  of 
bearings  are  small.  Whenever  these  angles  are  less  than  30  degrees,  the  bias 
behavior  is  such  that  if  the  angle  is  decreased  by  a  factor  of  n,  then  the 
bias  will  increase  by  a  factor  of  nl.  So,  if  the  angle  is  decreased  from  12 
degrees  to  6  degrees,  the  bias  will  increase  by  a  factor  of  four.  Bias  can 
also  be  a  problem  if  the  angles  between  the  true  lines  of  bearings  are 
sufficiently  close  to  1 80  degrees,  although  these  angles  will  not  generally 
cause  as  severe  a  problem  as  the  thin  angles. 

An  illustration  of  a  situation  in  which  the  magnitude  of  bias  is  large  is 
depicted  in  Figures  4  and  5.  In  the  first  diagram,  o  is  small,  and  yet  bias 

is  still  rather  large  due  to  the  thinness  of  the  angles  between  bearings.  The 

next  diagram  clearly  shows  the  Impact  of  a.  The  already  large  bias  has  been 
Increased  by  a  factor  of  nine,  making  it  enormous,  due  to  the  tripling  of  o. 

Note  that  increasing  o  by  a  factor  of  two  only  increases  the  lengths  of 
the  major  and  minor  axes  of  the  error  ellipse*  by  a  factor  of  two.  The  bias, 
on  the  other  hand,  Increases  by  a  factor  of  four.  At  some  point,  bias  will 
actually  reside  outside  the  ellipse.  The  relative  importance  of  random  error 
represented  by  both  the  ellipse  and  bias  changes  dramatically  with  changes  in 
c. 

In  the  general  cases  under  study,  however,  bias  appeared  to  be 
relatively  small.  As  long  as  the  angles  between  true  lines  of  bearings  are 

not  that  small,  then,  even  if  a  equals  3  degrees,  bias  will  not  be  very  large. 

See  Figure  6.  With  o«2  degrees,  the  "bias  region"  is  much,  much  smaller  than 
the  "spread  region"  ABCD  for  this  typical  configuration.  On  the  other  hand, 
if  the  angles  between  the  bearings  are  extremely  thin,  then  bias  can  be  large 
even  with  o«l  degree  (as  in  Figure  4). 

1  The  error  ellipse  is  a  confidence  ellipse  which  contains  the  emitter  with 
probability  (confidence)  8.  8  is  usually  chosen  to  be  .95 


As  mentioned  before,  the  general  expressions  for  first-order  bias  defy 
simple  qualitative  characterization.  In  their  general  form,  the  expressions 
contain  many  factors  which  contribute  to  bias  behavior.  Besides  the  two 
crucial  factors  previously  mentioned,  the  standard  deviation  and  the  angle 
between  LOBs,  the  next  most  significant  factor  which  contributes  to  bias  is 
sample  size.  So,  one  approach  of  analysis  has  been  to  examine  the  behavior  of 
bias  at  the  extreme  ends  of  sample  size:  the  2  LOB  case  and  the  infinite  LOB 
case. 

1 .  2  LOB  Case 

In  the  2  LOB  case,  all  of  the  fix  algorithms  will  produce  the  same 
estimate:  simply  the  intersection  of  the  2  LOBs.  Refer  to  Figure  7.  The 
direction  of  bias  is  determined  by  the  relative  location  of  the  emitter  to  the 
sensors.  If  the  emitter  is  outside  the  circle  shown,  then  bias  will  point 
away  from  the  sensors;  if  the  emitter  is  on  the  circle,  first-order  bias 
equals  zero;  and  if  the  emitter  lies  within  the  circle,  then  bias  will  point 
towards  the  sensors.  Generally,  however,  bias  has  been  analyzed  in  terms  of 
the  angle  between  the  true  bearings.  The  size  of  this  angle  is  directly 
correlated  to  the  emitter's  position  relative  to  the  circle.  When  the  emitter 
is  outside  the  circle,  the  angle  will  be  less  than  90  degrees;  inside  the 
circle,  greater  than  90  degrees;  and  on  the  circle,  the  angle  between  the  true 
bearings  will  equal  90  degrees. 

2.  Infinite  LOB  Case 

In  the  infinite  LOB  case,  bias  is  dependent  on  the  fix  algorithm  applied. 
Two  of  the  methods,  the  Minimization  of  Angular  Error  and  the  Sine  of  Error 
Minimization,  will  produce  no  first-order  bias  in  their  estimates  as  sample 
size  approaches  infinity  (asymptotically  first-order  unbiased),  while  the 
other  two  methods.  Perpendicular  and  Weighted  Perpendicular,  are  not  unbiased 
as  sample  size  approaches  infinity. 

It  is  worthy  of  noting  the  relative  Importance  of  bias,  as  opposed  to  its 
absolute  size,  in  the  2  LOB  case  as  compared  to  the  Infinite  LOB  case.  For 
the  Weighted  Perpendicular  and  the  Perpend icular  methods,  the  error  in  the 
estimate  for  the  2  LOB  case  (or  any  case  with  small  sample  size)  consists 
mostly  of  noise  and  bias  is  a  relatively  small  percentage  of  the  error.  As 
sample  size  approaches  infinity,  however,  the  noise  shrinks  (corresponding  to 
the  shrinking  of  confidence  ellipses),  and  the  error  in  the  estimate  is 
composed  almost  entirely  of  bias.  Also,  the  expected  location  of  the  infinite 
LOB  estimate  is  'short  of  the  true'  for  these  two  methods. 

For  the  Minimization  of  Angular  Error  and  the  Sine  of  Error  Minimization 
methods,  the  relative  Importance  of  bias  does  not  increase  as  sample  size 
increases.  Rather,  the  relative  importance  decreases  because  the  bias  shrinks 
faster  than  the  confidence  ellipse  as  LOBs  are  added. 

B.  Intermediate  Cases 

When  the  number  of  LOBs  is  greater  than  2,  the  "expected  location"  and 
hence  bias  lies  between  the  2  LOB  and  the  infinite  LOB  cases  (between  can  mean 
zero  bias,  however).  At  both  the  2  LOB  and  infinite  LOB  cases,  bias  has 
already  been  well-characterized.  In  between  these  boundaries,  the  important 


a 

aspect  to  note  about  the  properties  of  bias  is  the  movement  of  bias  between 
the  two  extreme  cases.  This  movement  betwee  the  extremes  is  also 
well-characterized  if  the  configuration  of  the  sensors  is  maintained  and  if 
more  readings  are  simply  taken  proportionately  from  each  sensor.  In  general, 
however,  in  order  to  characterize  bias  in  between  these  boundaries,  it  was 
necessary  to  restrict  attention  to  particular  configurations  of  emitters  and 
sensors.  These  particular  results  for  the  special  configurations  are 
discussed  in  the  next  section. 


* 
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III.  PROPERTIES  OF  THE  GENERAL  FORMULAS 


In  this  section,  various  results  regarding  bias  are  listed  according  to 
particular  cases.  Mathematical  formulas  ,  derivations,  and  proofs  are 
contained  in  the  Mathematical  Appendix. 

(Many  of  the  results  cited  in  this  section  have  been  described  under  the 
assumption  that  the  sensors  are  only  located  on  one  side  of  the  emitter) 

Four  general  2-dimensional  fix  algorithms  were  analyzed  for  the  purpose 
of  this  report.  Three  of  these  algorithms  (all  but  the  Weighted 
Perpendicular)  can  be  defined  as  Least  Squares  Fix  Algorithms,  which  are 
algorithms  that  estimate  parameters  by  minimizing  a  sum  of  squared  terms 
(usually  a  measure  of  error).  A  general  formula  for  least  squares  algorithms 
has  been  derived  (located  in  the  appendix)  which  is  applicable  to  all  three 
methods.  The  particular  formulas  then  only  differ  by  the  particular  term  that 
is  to  be  squared  and  summed.  For  example,  the  Perpendicular  Method  uses  the 
perpendicular  distance  from  the  estimate  to  the  LOB  as  the  term  to  be  squared 
while  the  Minimization  of  Angular  Error  Method  uses  the  angle  between  the  LOB 
and  the  line  from  the  sensor  to  the  location  estimate  as  the  term  to  be 
squared.  Since  the  formula  applies  to  all  three  least  squares  algorithms,  any 
analysis  of  the  formula  will  be  relevant  to  all  three  methods.  For  example, 
the  infinite  LOB  case  has  been  simply  characterized  and  also  the  movement  of 
bias  as  sample  size  increases  from  2  to  infinity  by  using  the  general  formula. 

The  results  regarding  bias  in  this  section  of  the  report  only  refer  to 
first-order  bias  and  not  to  total  bias.  When  the  standard  deviation  is  small, 
however,  the  higher  orders  of  the  bias  term  do  not  contribute  very  much  and 
first-order  bias  is  a  close  approximation  of  total  bias.  Simulation  results 
have  been  used  to  verify  that  this  approximation  is  reasonable. 

In  the  rest  of  this  section,  the  properties  of  each  of  the  algorithms 
except  for  the  Weighted  Perpendicular  method  are  described  in  depth.  For  the 
Weighted  Perpendicular  method,  there  is  an  extra  complication.  Relative 
distance  from  the  sensors  to  the  emitter  has  an  effect.  Except  for  this 
factor,  the  Weighted  Perpendicular  method  is  identical  to  the  Perpendicular 
method.  In  general,  the  Weighted  Perpendicular  method  emphasizes  the  data 
from  the  closest  sensors  while  the  Perpendicular  method  averages  the  data  from 
all  the  sensors.  Thus,  conceptually,  the  Weighted  Perpendicular  behaves  like 
the  Perpendicular  "for  the  closest  sensors". 

First-order  bias  for  the  Sine  of  Error  Minimization  method  is  identical  to 
that  given  by  the  Minimization  of  Angular  Error  method.  .So,  properties 
regarding  first-order  bias  for  the  Minimization  of  Angular  Error  method  can  be 
assumed  also  to  hold  for  the  Sine  of  Error  Minimization  method. 


A.  Extreme  Cases  --  Sample  Size 


1 .  2  LOB  Case 


In  the  case  when  there  are  only  2  LOBs,  the  four  fix  algorithms  will 
produce  identical  estimates.  The  estimate  is  simply  the  intersection  of  the  2 
LOBs.  The  general  scenario  of  the  2  LOB  case  is  depicted  in  Figure  8.  The 
true  bearings  are  shown  going  through  the  emitter  with  corresponding  angles  9, 
and  8j  while  e,  and  e2  are  the  angular  errors  associated  with  each  sensor. 

With  the  natural  assumption  of  symmetry  for  angular  error  (errors  of  e  and  -e 
are  equally  likely  —  this  assumption  is  used  throughout  this  report),  the 
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region  ABCD  may  be  constructed  as  shown.  LOBs  taken  in  this  situation  with 
angular  errors  smaller  than  the  bounds  e,  and  c2  will  Intersect  within  this 
region  and  consequently,  ABCD  is  a  geometric  representation  of  the  spread  of 
all  possible  estimates  (with  angular  errors  smaller  than  the  bound). 

Bias,  however,  is  much  smaller.  Using  the  previous  assumption  of 
symmetry  for  angular  error,  the  "region  of  bias"  formed  by  cancellation  of 
symmetric  error  in  ABCD  is  shown  back  in  Figure  6.  Furthermore,  it  is 
possible  to  show  that  this  "bias  region"  resulting  from  cancellation  of 
symmetric  error  will  always  be  wholly  situated  on  one  side  of  the  emitter, 
depending  on  the  size  of  the  angle  between  the  true  bearings.  Refer  back  to 
Figure  3*  The  angle  between  the  true  bearings  is  less  than  90  degrees  and 
thus  the  bias  region  is  located  entirely  on  the  far  side  of  the  emitter.  If 
the  angle  is  greater  than  90  degrees,  as  in  Figure  9,  then  the  bias  region 
will  be  entirely  on  the  close  side  of  the  emitter. 

2.  Infinite  LOB  Case 

In  the  case  of  an  Infinite  number  of  LOBs,  three  of  the  methods  can  be 
easily  analyzed  using  the  least  squares  general  formula.  This  analysis, 
however,  assumes  that  the  configuration  of  the  sensors  is  not  changed 
significantly  as  sample  size  goes  to  infinity.  If  LOBs  are  added  from  sensors 
so  that  the  total  geometry  of  the  situation  is  indeed  changed  significantly, 
then  the  behavior  of  bias  is  not  as  easily  predicted  (remember,  the 
configuration  can  have  a  major  effect  on  both  the  direction  and  magnitude  of 
bias).  The  analysis  for  the  least  squares  methods  is  as  follows. 

Least  Squares  Methods 

As  sample  size  Increases,  first-order  bias  for  the  Minimization  of 
Angular  Error  (and  Sine  of  Error  Minimization  )  method  approaches  zero.  In 
Figures  10-11 ,  a  comparison  is  made  between  a  particular  configuration  using  2 
LOBs  and  then  10  LOBs.  As  the  diagram  illustrates,  the  increase  in  sample 
size  caused  bias  to  shrink  (also  the  spread  of  the  fixes  decreased).  The 
speed  at  which  the  bias  will  shrink  can  be  affected  by  significant  changes  in 
the  configuration  as  LOBs  are  added  (including  changes  in  proportion  of  LOBs 
coming  from  a  certain  location).  If,  as  in  the  diagram,  there  are  no  changes 
in  the  configuration,  then  the  amount  by  which  bias  has  shrunk  is  proportional 
to  1/n  (n  being  the  factor  by  which  the  sample  size  increased). 

In  general,  asymptotic  first-order  bias  for  the  Perpendicular  method  can 
be  characterized  in  terms  of  the  location  of  the  emitter  relative  to  the 
average  sensor  location.  Refer  to  Figure  12.  The  heavy  dotted  line  through 
the  emitter  is  perpendicular  to  the  thin  dotted  line  from  the  average  sensor 
location  to  the  emitter.  As  sample  size  approaches  infinity,  the  expected 
value  of  the  estimated  location  of  the  emitter  will  lie  short  of  the  heavy 
dotted  line.  Without  oversimplifying  too  much,  one  may  summarize  by  saying 
that  for  sufficiently  large  sample  size,  one  expects  to  be  short  of  the  true. 

Weighted  Perpendicular  Method 

The  behavior  of  bias  for  the  Weighted  Perpendicular  method  from  the  2  LOB 
to  infinite  LOB  case  need  not  be  the  same  as  the  behavior  of  the  other 
methods.  Since  it  is  not  one  of  the  least  squares  methods,  it  cannot  be 
analyzed  in  the  same  manner  as  the  others.  However,  in  general,  asymptotic 
first-order  bias  for  the  Weighted  Perpendicular  method  will  behave  like  that 
of  the  Perpendicular  method. 

In  the  special  situation  when  there  are  only  two  symmetric  sensors, 
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however,  first-order  bias  is  identical  for  the  Perpendicular  and  the  Weighted 
Perpendicular  methods  (this  is  explained  on  the  next  page,  under  "Symmetric 
Case")  and  so  its  behavior  as  sample  size  increases  is  also  identical. 

B.  Transition  —  Sample  Size 

The  movement  of  bias  in  between  the  boundary  cases  for  the  three  least 
squares  methods  can  be  determined  by  examining  the  general  least  squares 
formula.  As  the  number  of  LOBs  Increases  from  2  to  infinity,  the  expected 
value  of  the  estimated  location  will  lie  along  the  straight  line  from  the  2 
LOB  estimate  to  the  infinite  LOB  estimate  (if  there  are  no  significant  changes 
in  the  configuration).  The  amount  by  which  the  expected  location  differs  from 
the  infinite  LOB  estimate  is  proportional  to  1 /n  (n  being  the  factor  by  which 
the  sample  size  increased)  See  Figure  13*  When  there  are  2  LOBs,  the 
expected  value  of  the  estimate  is  at  Point  A,  and  when  there  are  an  infinite 
number  of  LOBs,  the  expected  value  of  the  estimate  is  at  Point  B.  As 
illustrated,  the  convergence  is  very  fast;  at  only  4  LOBs,  the  expected  value 
of  the  estimate  is  half  the  distance  towards  the  infinite  LOB  estimate.  If 
sample  size  is  doubled  again,  the  distance  will  be  cut  in  half  again.  At  only 
16  LOBs,  almost  90X  of  the  distance  to  the  infinite  LOB  estimate  is  already 
covered. 

This  same  behavior  will  hold  even  if  the  starting  sample  size  is  not  2 
LOBs.  For  instance,  if  LOBs  are  initially  added  to  the  first  2  LOBs  such  that 
the  configuration  is  changed  significantly,  then  the  above  analysis  will  not 
be  valid.  If  the  configuration  does  become  stable  at  some  point,  so  that  all 
additional  LOBs  will  not  seriously  affect  the  configuration,  then  the  bias 
will  begin  to  converge  from  that  point  towards  the  infinite  LOB  estimate.  For 
example,  assume  that  the  first  4  LOBs  have  been  attained  such  that  each  LOB 
changed  the  configuration  and  thus  bias.  Now,  each  additional  LOB  after  the 
4th  LOB  will  be  added  without  significant  changes  in  the  arrangement.  Then, 
the  expected  value  of  the  estimated  location  will  now  move  in  a  straight  line 
towards  the  infinite  LOB  estimate,  moving  half  the  remaining  distance  with 
each  doubling  of  sample  size.  Thus,  at  8  LOBs,  the  expected  value  of  the 
estimate  will  be  half  the  distance  towards  the  infinite  LOB  estimate. 


C.  Variation  of  Bias  With  Position  —  Special  Cases 

The  following  special  cases  have  been  included  in  order  to  provide 
certain  Insights  to  the  behavior  of  bias. 

1 .  Symmetric  Case 

The  general  symmetric  scenario  exists  when  the  sensor  configuration  is 
symmetric  about  a  line  through  the  emitter;  that  is,  the  configuration  on  one 
side  is  a  mirror  image  of  that  on  the  other.  In  addition,  corresponding 
sensors  on  either  side  have  been  assumed  to  be  equal  in  accuracy. 

If  it  is  assumed  that  corresponding  sensors  are  equal  in  accuracy,  then 
bias  will  lie  entirely  along  the  line  of  symmetry,  regardless  of  the  method. 
Thus,  the  symmetric  assumption  was  chosen  as  the  natural  method  of  reducing 
bias  to  the  determination  of  extent  short  and  long. 

Minimization  of  Angular  Error  Method 

For  the  Minimization  of  Angular  Error  method,  first-order  bias  in  the 
symmetric  situation  is  easily  characterized.  See  Figures  14-15.  If  ail  of 


the  angles  between  the  true  bearings  from  symmetric  sensors  are  less  than  90 
degrees,  then  bias  will  point  long  (along  the  line  of  symmetry)  and  if  all  of 
the  angles  are  greater  than  90  degrees  and  less  than  1 80  degrees,  then  bias 
will  point  short.  And  if  the  angle  equals  90  degrees,  then  first-order  bias 
will  equal  zero. 

Perpendicular  Method 

In  a  symmetric  situation,  Perpendicular  bias  has  been  defined  with 
respect  to  the  angle  between  true  bearings  and  sample  size.  See  Figures 
16-19.  When  there  are  only  2  LOBs,  the  bias  will  point  long  along  the  line  of 
symmetry  when  the  angle  between  the  true  bearings  is  less  than  90  degrees;  it 
will  point  short  when  the  angle  is  greater  than  90  degrees;  and  the  bias  will 
be  zero  when  the  angle  between  the  bearings  equals  90  degrees  (this  behavior 
is  consistent  with  that  already  shown  for  the  general  2  LOB  case). 

As  LOBs  are  added,  both  the  direction  and  the  magnitude  of  the  bias 
vector  can  change  depending  on  the  size  of  the  angle  between  the  true  lines  of 
bearing.  In  Figures  16  and  17,  all  of  these  angles  are  greater  than  120 
degrees.  Bias  is  already  short  using  this  configuration  when  there  are  only  2 
LOBs;  as  sample  size  increases,  bias  will  shrink  and  become  less  short  (along 
the  line  of  symmetry). 

If  all  of  the  angles  between  true  bearings  are  equal  to  120  degrees,  then 
the  bias  will  not  change  in  size  or  direction  as  sample  size  Increases. 

If  these  angles  are  less  than  120  degrees  (as  in  Figures  18  and  19),  the 
bias  vector  will  point  short  from  whatever  estimate  it  pointed  to  when  there 
were  only  2  LOBs.  So,  if  the  bias  pointed  long  at  2  LOBs  (angles  less  then  90 
degrees),  then  it  will  change  direction  and  point  short.  And  if  the  bias 
pointed  short  at  2  LOBs  (angles  greater  than  90  degrees  but  less  then  120 
degrees),  then  bias  will  point  even  shorter  and  increase  in  magnitude. 


Weighted  Perpendicular  Method 


When  the  relative  distances  from  the  sensors  to  the  expected  value  of  the 
estimated  location  are  all  the  same,  then  the  Weighted  Perpendicular  method 
will  produce  exactly  the  same  first-order  bias  as  the  Perpendicular  method. 
Thus,  when  all  of  the  LOBs  come  from  only  two  symmetric  sensors,  first-order 
bias  will  be  identical  for  both  of  the  methods. 


2.  Effects  of  Range  Differences  Between  Sensors 


How  each  method  uses  the  data  it  receives  has  also  been  studied  to  some 
extent. 


Perpendicular  and  Weighted  Perpendicular  Methods 


As  an  example  of  how  the  Perpendicular  and  the  Weighted  Perpendicular 
methods  can  deviate  in  certain  cases,  a  special  configuration  was  used.  In 
Figures  20-21 ,  the  sensors  are  placed  on  two  concentric  arcs  so  that  sensors 
on  the  same  arc  are  equidistant  from  the  emitter.  Furthermore,  a  sensor  on 
one  arc  can  be  found  placed  exactly  in  line  with  the  emitter  and  with  a  sensor 
on  the  other  arc  (so  that  the  angle  between  the  emitter  and  consecutive 
sensors  are  the  same  from  one  arc  to  the  other).  Thus,  the  only  difference 
between  the  sensors  on  one  arc  from  those  on  the  other  is  the  distance  from 
the  emitter. 

As  the  figures  illustrate,  there  is  a  marked  difference  between  the 


estimate  that  each  method  produces.  The  Weighted  Perpendicular  method  has 
used  the  data  from  the  closest  sensors  (and  ignored  the  data  from  the  farthest 
sensors)  and  has  come  up  with  a  better  estimate  while  the  Perpendicular  method 
has  averaged  the  data  from  the  sensors  and  has  come  up  with  a  worse  estimate. 

Minimization  of  Angular  Error 

In  an  Identical  configuration  as  the  one  described  above,  the 
Minimization  of  Angular  Error  method  will  weight  the  data  from  the  two 
concentric  arcs  so  that  only  the  data  from  the  closest  sensors  is  used  and  the 
data  from  the  farthest  sensors  is  virtually  ignored. 

3.  Effects  of  Angles 

Minimization  of  Angular  Error 

The  general  90  degree  case  exists  when  there  are  only  two  sensor 
locations  and  the  emitter  is  located  such  that  the  true  bearings  are  90 
degrees  apart  (See  Figure  22).  This  special  case  is  a  generalization  of  the  2 
LOB,  90  degree  case  (illustrated  in  Figure  7),  for  which  first-order  bias 
equals  zero.  When  multiple  LOBs  are  taken  from  the  sensors,  first-order  bias 
using  the  Minimization  of  Angular  Error  method  equals  zero,  regardless  of  the 
respective  number  of  L0B3  taken  from  each  sensor  location  or  of  the  angular 
standard  deviation  associated  with  each  sensor. 

Perpendicular  and  Weighted  Perpendicular  Methods 


The  Perpendicular  and  Weighted  Perpendicular  methods  are  not  unbiased 
when  the  true  bearings  are  90  degrees  apart  except  when  there  are  only  2  LOBs. 
For  more  than  2  LOBs,  the  bias  for  the  Perpendicular  method  Is  not  easily 
characterized  because,  unlike  the  Minimization  of  Angular  Error  method,  its 
bia3  depends  on  the  respective  number  of  LOBs  from  each  sensor  and  the  angular 
standard  deviation.  As  sample  size -approaches  infinity,  however,  the  estimate 
produced  by  either  method  will  be  biased  short  of  the  true  location. 
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FIGURE  7 
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FIGURE  10 
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FIGURE  11 
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FIGURE  12 
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Expected  Value  for  the  Estimated  Location 

(movement  with  respect  to  sample  size) 


This  analysis  is  valid  for  three  of  the  methods  (all  but  the  Weighted 
Perpendicular).  Also,  the  result  assumes  that  as  N  increases,  there 
are  no  significant  changes  to  configuration. 
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FIGURE  16 
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FIGURE  17 
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MATHEMATICAL  APPENDIX 


Formulas  for  bias  have  been  obtained  for  the  four  algorithms  discussed 
in  this  report.  Three  of  the  algorithms  have  been  characterized  by  least 
squares  while  the  bias  formula  for  the  Weighted  Perpendicular  method  has  been 
derived  separately  at  the  end  of  this  section. 


2-Dimensional  Least  Sauares  Fix  Algorithms 


Least  Squares  Algorithms  are  algorithms  that  estimate  parameters  by 
minimizing  a  sum  of  squared  terms  (usually  a  measure  of  error)  where  each  term 
depends  only  on  the  parameters  and  one  observation.  For  our  application  the 
parameters  are  the  x  and  y  coordinates  of  the  location  estimate  and  the 
observations  are  LOBs. 


Definitions 


( X, Y)  -  true  location  of  the  emitter 


0^ ( X, Y)  -  true  bearing  from  the  kth  sensor  to  the  emitter  (X,Y) 


0,  -  observed  bearing  from  the  kth  sensor  (multiple  readings 


are  treated  as  coming  from  different  sensors) 


(x,y)  -  (x(0,...,0  ))«estimated  location  of  the 

1  n  1  n  emitter 


standard  deviation  of  the  angular  measurement  of  the 


kth  LOB  in  radians  (multiply  by  ir/180  if  in  degrees) 


ek  -  error  in  kth  LOB  ( ©^-©^(X, Y) ) 


(xk,yk)“  sensor  location  of  kth  LOB 


rl  -  (x-x.  )2+(y-yJZ  R?-(X-x.  )2+(Y-y.  )2 


*  "V  '•  /(  '«  .  •  ’I  '»  m~*.  /■  ,  •/  *  .*»'**'*  •  h  "»  *  ", 


>  >  *  ■  /  V  .* 


r 


LIKELIHOOD  MINIMIZATION  TYPE  FIX  METHODS 


Part  of  the  analysis  of  these  methods  is  independent  of  which  method  is 
being  analyzed. 

Extra  definitions  that  apply  to  Likelihood  based  methods 

L^  -  Lk(xfy,9k)-  the  'squared  error  term'  corresponding  to 

the  kth  LOB  and  a  location  estimate  of  (x,y) 
for  the  least  square  method. 

A  A 

LCx.y.e^ . 9^)  -  ^I^L^-the  sum  of  squares  that  the  method  minimizes 

Assume  that  x,y  are  defined  implicitly  as  the  minimum  of 

A 

L-E  Lk(x,y,8k)  (for  summation  notation,  see  *  below) 

in  the  sense  that  3L/3x-dL/3y-0 

The  evaluated  partial  derivatives  of  L^  will  be  denoted  as  follows 
(recall  that  evaluated  partial  derivatives  are  used  in  a  Taylor  Series). 

Let  bk»<3Lk/dx  ck»3Lk/3y  dk«32Lk/3x2  e^-d2L^/dxdy  fk«d2Lk/dy2 

gk-33Lj</3x3  hk»33Lk/3x2  y  ik-33Lk/c9x3y2  j^^L^y3 

Let  '  notation  represent  derivatives  with  respect  to  ek  for  bk  to  fk> 
Then,  since  3L/3x-0,  Zb  -0  and  since  3L/dy-0,  Ec  -0 

K  K 

O-3/3e1(3L/3x)-E3/3ei*bk-EOx/3e1*3bk/3x+3y/<3ei*3bk/<3y+30k/3Ei*i3bk/39k) 

-  3x/3ei(Edk)+3y/3e1(Eek)+b^ 

0-32/3e2(3L/3x)-  32x/3e2(Edk)*(f?x/3ei)  2(E3dk/3x)  +  ((3x/3ei)  ( 3y/3ei)  (E3dk/3y) 
+3x/3e1*d|+32y/3cJ(Eek)  +  Ox/3e1)  Oy/3e1)(E3ek/3x)  +  (dy/3e1)  2(Edek/<3y ) 
+3y/3e1»e|-*-dx/3e1*8b^/3x+3y/3e1(3b^/3y)+b^ 
-32x/3e2(Edk)+(3x/3ei)2(Egk)^(3x/3ei)(3y/3e1)(Ehk)+3x/3e1*d^+32y/3e2(Eek) 
♦(3x/3ei)(3y/3ei)(Ehk)  +  (3y/3ci)  2(Elk)+3y/3ci*e^3x/3ei*3b^/3x 
♦3y/3£i(3b|/3y)*b^ 

*  n 

Un-labeled  sums  are  assumed  to  be  (n  equals  the  number  of  LOBs). 

Sums  labeled  E^  are  assumed  to  be  ^ g” 


•d2x/de2(  Ed  )  +32y/3e?  ( Ee  )  +  Ox/3e  )  2(  Eg.  )  +20x/<3e  . )  Oy/c3e  . )  (Eh  ) 

IK  IK  IK  i  ik 

>(ay/ae.)J(Eik)^2(ay/aei)e^20x/aE.)d!*b^ 


Similar  argument  may  be  applied  to  3L/3y  yielding 


d/dc . 


3L/3x 


1  3L/3y 


Edk  Eek 
Eek  Efk 


dx/di 


3y/3ei 


f 

i  bi 
♦ 

» 


d  y/3  e 


( Ed.  )  ( Ef,.)-( Ee,  )  2  -Ee, 


3L/3x 

3L/3y 


Edk  Eek 
Eek  Efk 


a*x/ae2 

32y/3e2 


TERM  1 


TERM  2 


where 

TERM  ll  r(ax/aei)2(Egk)+2(ax/aci)(ay/ae1)(Ehk)+(ay/ae1)I(Elk) 


TERM  2j  [(ax/aei)2(Ehk)+2(ax/aei)(ay/ae1)(Eik)+(3y/a£1)2(Ejk) 


2(3y/3ei )e j+2(3x/3ei )dj+bj 
2(ay/3e1)f j+2(3x/aE1)ej+c^ 


Thus, 


a2x/3e2 

a2y/3e2 


(Edk)(Efk)-(Eek)2 


-Eek  TERM  1 


Edk  TERM  2 


31x/3e2 

First-order  bias-  1/2  •  Z  -a2 

1  aay /del  1 


Particular  algorithms  use  a  different  Lk  from  which  to  derive  the  partial 
derivatives  used  in  the  formula  for  bias. 

Perpendicular  Method 


Lk-  [ (x-xR) 2cos28k+(y-yk) 2sin2ek-2(y-yk) (x-xk)sinekcos0k] 


The  evaluated  partial  derivatives  are  as  follows  (  evaluated  at  8  »e  ) 

k  k 

b’-  -2( Y-y±) 

c!-  2(X-xi) 

d’-  -U(Y-yi)(X-x1)/R* 

ej-  2((X-x. )2-(Y-y.)2)/R2 

f’-  ^ ( Y-y . ) (X-x  ) /R2 
1  i  i  l 

b%  4(X-x.) 
c%  *4 (  Y-y  ^ ) 

V  2U-yk)VR' 

•k-  -2(v-ykxx-*k>/Rj 


rk- 

8k  ■  \  •  ‘k  ■  Jk  -° 

Minimization  of  Angular  Error  Method 

L^-  [Arctani (x-xk)/(y-yk) )-0k]2 

The  evaluated  partial  derivatives  are  as  follows  (  evaluated  at  ek-9k) 

b|-  -2( Y-y t ) /R2 

c|-  2(X-Xi)/R2 

d[-  4(Y-yi)(X-xi)/Rj 

e’-  2((Y-y1)2-(X-xi)2)/Rj 

fj-  -4(Y-yi)(X-x1)/Rj 

It 

v° 

ft 

ci-  0 

V  2(l-V2/R!! 

V  -2(l-rkHX-,t)Zi; 

fk-  2<«-V2/< 

V  -'2<x-«k)(x-yk)'2sk 

hk-  (8(X-«k),(X-yk)-»(Y-yk)3)/RR 
1,  -  ( 8 ( X-x,. )  ( Y-y,. )  2-4  ( X-x.  ) 3 )  /R^ 


*1. 


^Xwjr*JT'Jrmjr*Ji  V  ’/•  ^  V  V.V  /^.  r.  tr j  "• 


jk-  -12(X-xk)i(Y-yk)/R6 
Sine  of  Error  Minimization  Method 

The  partial  derivatives  for  this  method  are  identical  to  the  ones  just 
shown  for  the  Minimization  of  Angular  Error  method.  Thus,  first-order  bias  is 
identical  for  both  methods.  However,  the  Lk  is  as  follows. 

Lk«  [ (x-xk) Icos2ek  +  (y-yk) 2sin20k~2 ( y-yk ) (x-xk) sin0kcos0k ]/[ ( x-xk ) 2  +  (y-yk) 2  ]  ^ 

Notice  the  similarity  of  this  term  with  the  Lk  for  the  Perpendicular 
Method.  It  is  the  same  except  that  this  one  is  divided  by  (x-xk)  2+(y-yk) 2 , 
which  is  simply  r2,  the  distance  squared  from  the  kth  sensor  to  the 
estimated  location  of  the  emitter. 

Infinite  Number  of  LOBs 

First-order  bias  as  the  number  of  LOBs  approaches  infinity  is  easily 
characterized.  Each  of  the  terms  in  the  general  formula  has  a  sum  in  the 
numerator  and  in  the  denominator,  except  for  the  last  terms  in  Term  1  and 

h  n 

Term  2:  b^  and  c^.  As  N  goes  to  infinity,  those  terms  with  sums  in  both 

numerator  and  denominator  will  go  to  zero  and  the  asymptotic  first-order  bias 
will  be: 


X 

-Ee 

;Iek 

Ed 

Udk)(lfk)-(rek)2  |_~lek  ldkJ  |^J 

For  the  two  methods  Minimization  of  Angular  Error  and  the  Sine  Of  Error 

Minimization,  the  b^  and  the  c^  terms  are  zero  so  the  two  methods  are 

asymptotically  unbiased  for  first-order. 

The  Perpendicular  method,  however,  is  asymptotically  biased.  When  the 
proper  terms  are  substituted  into  the  formula,  the  result  is  as  follows. 

With  a  change  of  coordinates,  the  y-axls  can  be  made  to  go  through  the 

emitter  and  the  average  sensor  location,  and  b^  can  be  made  to  equal  zero 

(without  loss  of  generality).  Then,  the  x-component  of  bias  will  be  ignored, 
and  only  the  y-component  examined.  The  y-component  of  asymptotic  first-order 
bias  will  be  as  follows: 


1/2*  E. 


(Edk)(ifk)-(Eek)2 


Ed  •  c  •  a2 
k  i  °i 


The  sign  of  this  term  will  always  be  negative  so  first-order  bias  will 
point  towards  the  average  sensor  location. 

A  similar  analysis  may  be  performed  in  order  to  determine  asymptotic 
first-order  bias  for  the  Weighted  Perpendicular  Method.  Like  the 
Perpendicular  Method,  the  Weighted  Perpendicular  Method  produces  biased 
estimates  as  sample  size  approaches  infinity.  In  general,  this  first-order 
asymptotic  bias  will  point  short  of  the  true  location  (details  are  discussed 
in  the  Weighted  Perpendicular  section). 
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General  Symmetric  Case 

Symmetry  of  sensors  produces  the  following 
Eek  •  ■  e1k  •  0 

Then,  <3x/3e^  -  -bj/Ed^  and  3; 

and, 


Then,  dx/de^  -  -bj/Id^ 


3y/3ei  -  ~cVEfk 


32x/3e2 

■ 

i 

m 

32y/<3e2 

If  it  is  also  assumed  that  the  opposite  pairs  of  o^s  are  the  same,  then 


all  of  the  terms  In  the  top  component  of  the  second  matrix  will  be  zero  after 
they  are  summed  on  i.  Then,  the  x-component  of  bias  will  be  zero  and  the 
y-component  will  be  as  follows. 


Symmetric  bias  (y-component)-  1/2  •  Zi(32y /de*)  •  a2  where 
3*y/a£2  -  (-bi/Zdk)2(Zhk)-*-(-c./Efk)2(ZJ1<)+2(-ci/Zfk)f^2(-bi/Zdk)e^cJ 


If  less  assumptions  are  made  (for  instance,  If  the  opposite  pairs  of 
are  not  the  same),  then  intermediate  results  may  also  be  derived  which  are1 


useful. 


Perpendicular  Method 


It  will  be  assumed  here  and  for  the  rest  of  the  symmetric  cases  that 
opposite  pairs  of  o^s  are  the  same.  Then,  the  x-component  of  first-order 

bias  is  zero  and  the  y-component  of  first-order  bias  is  as  follows. 


Ei[a2y/3e2]-o2  -  Zi[-o2/Zfk(-2-b[-e[/j:dk  -2-f[*c[/Zfk  +cj  )] 


With  more  than  2  LOBs,  bias  is  always  negative.  However,  as  sample  size 
increases,  the  direction  that  bias  will  move  in  is  dependent  upon  the  size  of 
the  angles  between  the  true  lines  of  bearings.  Bias  will  become  less  negative 
if  the  angles  between  bearings  are  all  greater  than  120  degrees  and  bias  will 
become  more  negative  if  the  angles  are  less  than  120  degrees.  If  all  of  the 
angles  are  equal  to  120  degrees,  then  bias  will  not  move  from  its  2  LOB 
position  to  its  Infinite  LOB  position.  The  proof  of  this  involves  determining 
the  sign  of  the  bias  term  when  the  angles  between  bearings  are  either  less 
than  or  greater  than  120  degrees. 

When  the  absolute  value  of  all  of  the  angles  of  the  true  bearings  from 
each  sensor  are  greater  than  or  equal  to  120  degrees, 


t  *  -i 

biei"  “i,r^C039i~C03j9^ 


cjfj-  8ri(cos8  . -cos^e^ )  cos^S  1/2 


V  1/2 


V  3/2 
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Zd  S  N/2  £  1/3*  3N/2  S  1/3‘E^ 

£fk>  3Zdk  so  1/Zfk$  1/(3*Edk)  and  -1/Zfk*  -1/(3*^) 

-*;'£v=;^v  -3v;/<3-cv  -  =;r;/(3-£v 

the  right  side  of  inequality  -(iJcose^  -IbcosV  )/(3»Zdk)  2  0 

Bias  has  one  more  negative  sign  in  it,  therefore  the  above  term  is  less 
than  or  equal  to  zero. 

In  order  to  find  out  the  movement  of  bias  as  sample  size  increases,  for 
the  opposite  case  when  all  of  the  angles  are  less  than  or  equal  to  120 
degrees,  the  same  argument  is  applied  except  that 

cose^  1/2  dk*  1/2  fkS  3/2  and  -1/Zfk*  -1/(3*Edk> 

so  that  bias  is  found  to  be  greater  than  or  equal  to  zero.  Combining  these 
two  results  provides  for  the  conclusion  of  its  behavior  at  120  degrees;  the 
movement  is  zero. 


Minimization  of  Angular  Error  Method 

After  insertion  of  the  proper  terms,  first-order  symmetric  bias  for 
sensors  with  the  same  standard  deviation  will  be  as  follows: 

2o2  ZZ  (X-xk)2(Y-y|<)[(Y-yj)2-(X-x1)2]/(Rk-R^) 

Z(X-xk)2/Rj  ZZ  (Y-yk)4(X-xJ)*/(Rj*R'j) 

The  sign  of  this  expression  depends  on  if  the  angles  between  bearings 
are  less  than  or  greater  than  90  degrees.  At  90  degrees,  first-order 
symmetric  bias  is  zero. 


Weighted  Perpendicular  (Not  Definable  in  Terms  of  Least  Squares) 

The  Weighted  Perpendicular  Method  is  not  definable  as  a  least  square's 
method  and  so  a  separate  formula  needs  to  be  derived  for  first-order  bias.  I 
is  necessary  to  establish  notation  for  this  case. 


Definitions 


-  cos20  ^  -  cosl8i(X, Y)-( Y-y^) 2/R* 

b  -  sin*8.  B.  -  sln20 . (X, Y)-(X-x  ) */R* 


i 


i 

M  A 

sine^ose^ 


i 


Ci-(X-xi)(Y-y1)/R* 


di  -  Vi  -  Vi 
ei  ■  -Vi  +  Vi 


°i  •  Vi  •  ciyi 

Ei  ■  -C1X1  *  Vi 


£vri 


a 

b 

c 

d 

e 


Eai/rI 

EVr’i 

Eci/rI 

Zdi/r2 

Iei/ri 


A  -  ZA./R2 
i  l 


B 

C 

D 


ZB./R2 

l  l 
ZC^R2 

ZD. /R2 

l  l 


E  -  ZE^R2 


•  - 

a  -c 

-c  b 

* 


* 


(x.y)  -  (Is./r 2)  '(Es./r2 


First  Partials: 

jr*/<30k  *  2(x-x1)(<3x/^ek)  ♦  2(y-y.)(<3y/30k) 

3a/3ek  -  -2ck/r*  -  Z(ai/r^*  dr*/dQk) 

3b/3ek  -  2ck/rk  ~  I(bi/ri*  ^ri/^®k) 

3c/3ek  -  (ak-bR)/r*  -  E(c./rJ.  3r2/30k> 

3d/3ek  -  C-2ckxk  ♦  <  Vak)yk]/rk  '  z(  Vri*  <^ri/a®k) 

3e/30k  -  [(bk-ak)xk-  2ckyk]/r*  -  ZCd^rJ-  3r*/3ek> 
Differentiating  the  defining  equation: 


3x/3o, 

I* 

dy/de. 


(a(3b/30,)-*-b(3a/30  >-20(30/30,.) ) 

-  r2\2  K  K  K 


(ab  -  c^) 


bd  +  ce 


cd  ♦  ae 


b(3d/3e,  )  ♦  d(3b/3e,  )  ♦  c(3e/3e  )  ♦  e(3c/30.  ) 


(ab  -  c2)  j^c(3d/30k>  +  d(3c/30k>  ♦  a (de/dQ^)  *  e(3a/30k> 

Simplifying  the  defining  equation: 


Let  fk«  2[b2d+bce-ace-c2d]ck  ♦  [2cbd+c2e+abe](ak-bk) 

♦  [b2a-c2b](-2ckxk  ♦  (bk'ak)yk)  *  !>(ab-c2) (Vak)xk  +  2ckyk) 

h  -  [-b2d-bce]a  ♦C~ace-c2d]b  ♦C2cbd*c2e^abe]ck+[b2a-c2b]dk+Cc(ab-c2) ]e! 

k  k  k 

Let  gk«  2Cbcd>c2e-acd-a2e]ck  ♦  [2ace-*-c2d>abd](ak-bk> 

♦  [c(ab-c2> ](-2ckxk  ♦  (bk-ak>yk>  ♦  [a(ab-c2> ]( (bk~ak)xk  ♦  2ckyk^ 

i  ■  [-bcd-c2e]a.  +[-acd-a2e]bk+[2ace*c2d+abd]Ck*[ (ab-c2)c]dk+[a(ab-c  ) ]e 

k  k  ^ 

3x/3©k  -  {fk/r*  -  Z  1/rj-  3r*/3ek-  hj  }/(ab-c2)2 

3y/30k  -  (gk/r*  “  E  1/rj*  ^rj/^ek*  lj  }/(ab-c2)2 

Let  Fk,Gk,Hk,  and  Ik  be  fk,gk,!V  and  *k  evaluated  at  the  true  33  in  the 
Taylor  Series  expansion. 


Plugging  in  for  3r2/30k  and  evaluating  at  the  true  yields 

dx/3ek  -  (Fk/R*  -Jl"2{((X-xj)3x/<30k  +  (Y-yj)3y/a0k)HJ/Rj}/(AB-Ci)J 
<3y/<30k  -  {Gk/Rk  -j£"2{((X-xj)ax/aek  +  (Y-yj)ay/a0k)Ij/Rj}/(AB-C2)2 
where  the  partials  shown  are  also  evaluated  at  the  true. 


q  -  C(AB-C2)2+2JZIj,((X-xj)Hj/Rj}]C(AB-C2)2*2jr!|l{(Y-yj)IJ/R'j}] 

-C2 j  Elj1  { ( Y-yj  ) H j /R J }  ][ 2  { ( X-x  )  I  /rJ  }  ] 

For  partials  evaluated  at  the  true 

3x/9 0k  -  [  F^R2  { ( AB-C 2 )  2  ♦  j  Z?2  { ( ( Y-Yj ) I j /R J } -Gk/R2 { j ^2 ( Y-Yj ) H  /R J } ] /q 

3y/3ik  -  C-F|</R2{jz!j»2{((X-xj)IJ/R'j}+Gk/R2{(AB-C2)2+jr!|l2(X-xj)Hj/R^}]/q 

The  answer  must  be  invariant  under  translations  (Independent  of  where  the 
origin  is  placed.)  Fortunately  calculations  are  a  good  deal  easier  if  we 
choose  (X, Y)-(0, 0) .  (Note  we  don't  do  that  in  the  computerized  examples 
however.)  With  this  assumption 

Dk-Ek-0  for  all  k.  Hence  D-E-0.  Further  implying  Hk~Ik-0  for  all  k. 
Hence  q»(AB-C2)**  and 

ax/90k-(AB-C2)2Fk/(R2q)-{B(-2Ckxk>(Bk-Ak)yk)+C((Bk-Ak)xk^2Ckyk)}/[R2(AB-C2)] 


k  k  '  k  k  '  k 


k  k'  k  k  k  k' 

i /rB2f  *R_p2'\  i 


-  {B(-(Bk+Ak)y|<)+C((Bk>Ak)xk}/[R^(AB-C‘i)] 

-  {-Byk+Cxk)/[R2(AB-C2)] 

3y/3ek-(AB-C2)2Gk/(R2q).{C(-2Ckxk>(Bk-Ak)yk)^A((Bk-Ak)xk>2Ckyk)}/CR2(AB-C2)] 

-  {-Cyk>Axk}/[R2(AB-C2)] 


Second  Partials: 

0lr*/aek-  2(x-xi) (32x/a02)*2(y-yi)3*y/a02+2(9x/<30k) *+2(3y/a0k) 1 

32a/3e2  —  2(ak~bk)/r2  +  ( 4ck/rj|) (<3r2/30k)  ♦  EUa^r®)  (3r2/90k) 2 

■►E(-a1/r^)(a2rJ/802) 

32 b/90*  -  2(ak-bk)/r2  -  (4ck/rJ)(dr2/30k)  *  I^bj/rJ)  (drj/^) 2 
♦E(-b,  /r!*)  (<32r?/<30,2) 


<32c/3e*  “-4ck/r£  "  2(  (ak-bk)/r4)(3r2/30k)  ♦  Z(2c  .  /r^)  2 

+Z(-c./r})(32r2/3e2) 

a2d/3e2  -  (Vk'2(Vbk)xk)/rk  +  2((1-2bk)yk+2ckxk)/rk)(ark/a®k) 

♦  Z(2di/r^) (3r*/3ek) 2  +  EC-dVr^)  (32r?/302) 

32e/3e2  -  (iJckxk>2(ak-bk)y(<)/r2  +■  2( ( 1 -2bk)xk-2ckyk)/rJ) (3r2/30k> 
+  E(2e  /r^)(3r2/30,)2  ♦  E(-e . /r4) (32r?/302) 

11  Ik  1  1  1  K 


32x/3e2  i 

k  -  T1  (3a/3@.  )2*T2(3a/3S.  )(3b/3S.  )*T3(3a/30.  )(3c/30.  ) 

a*y/3e2  (ab-c2)3  k  k  k  k  k 

+T4(3a/30k)(dd/30k)+T5(3a/30k)(3e/30k)+T6(3b/30k)2*T7(db/30k)(3c/30k) 

+T8Ob/a0k)(ad/<30k)+T9(3b/30k)Oe/a0k)+Tio(ac/30k)2*Ti  i  (ac/3ek)(ad/a0k) 

+T12(3c/d0k)(3e/a0k)*T13(d2a/ae2)*Tl4(32b/a02)+T15(32c/302) 


+T1 6( 32d/3i2)+T1 7(32e/3@2) 


where  T1  through  T1 6  are  as  listed  below. 


(no  more  derivatives  will  be  taken  so  evaluated  derivatives  will  be 
used  in  place  of  unevaluated  ones.  In  other  words,  upper-case  letters  will 
replace  lower-case  letters.) 


T1 -2B(BD+CE) 
T2- 
T3- 

L~°' 

T4-  (AB-C2) 

15-  (AB-C2) 


f  i 

B 

_C_ 


2ABCE*4C2BD-*-2C2EC 

2ABCD+4ClAE+2C2DC 

-8B2CD-6C2BE-2AB 

-6C*BD-4ABCE-2AB 

-2B* 

-2BC 

-2BC 
-2C* 


^  J 

2E  1 

2D-4C2EC  J 


T6-  2a'[c§:»e]-2<#b-c‘>[aS] 

n-  -8ic[?S:Se]-;<ab-c‘>[#eadC0] 

T8-  (AB-C2) 


T9-  (AB-C2) 


-2C* 
-2  AC 

-2  AC 
-2A2 


T10-  2C  2  f38D+CE 
[2CD+3AEJ 


+2AB  rBD^3CE] 
2AB  [3CD+AEJ 


4CB 

2AB+C2 

2AB-*-C2t 
4CA  J 


Til-  (AB-C2) 
T12-  (AB-C2) 
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T13- 

(AB-C2) 

T1  4- 

(AB-C2) 

T15- 

(AB-C2) 

T16- 

(AB-C2) 

T1  7- 

(AB-C2) 

Again, 

,  if  it 

-B2D-BC£ 

-BCD-C2E 

-ACE-C2D 

-ACD-A2E 


►ABE  I 
►  ABDj 


•LaJ 


T1 -T2-T3-T6-T7-T1 0-T1  3-T1  4-T1  5-0 


then 

(because  D-E-0) 


Infinite  Number  of  LOBs 


As  in  the  other  methods  studied  in  this  report,  the  first-order  bias 
vector  for  the  Weighted  Perpendicular  method  points  along  a  line  segment  from 
the  bias  for  the  2  LOB  case  to  the  bias  for  the  infinite  LOB  case  (if  there 
are  no  significant  changes  in  the  configuration).  First-order  bias  as  sample 
size  approaches  infinity  is  as  follows. 


v 

V 


v. 


■v. 


Asymptotic  first-order  bias- 

B  C 

C  A 

Z 

V 

-yi- 

•  °Vr\ 

AB 

-  C 

2 

V 


It  is  possible  to  set  C-0  by  finding  the  "average  angle"  of  the  angles 
between  LOBs  and  true  north.  Along  this  "average  angle",  asymptotic 
first-order  bias  will  lie  short. 


‘  ^  •-* ■" \ ft  m  A  k  fc. kV k  k.\  k  *  k^Vj 


